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Abstract: We consider the following trace function on n-tuples of positive operators: 

%(A ll A 2l ... 1 A n )=Tr((J2 AP 3 ) 1/P ) 

\ j= i / 

and prove that it is jointly concave for < p < 1 and convex for p = 2. We then derive from 

this a Minkowski type inequality for operators on a tensor product of three Hilbert spaces, 

and show how this implies the strong subadditivity of quantum mechanical entropy. For 

p > 2, $ p is neither convex nor concave. We conjecture that $ p is convex for 1 < p < 2, 

but our methods do not show this. 
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I. INTRODUCTION 

Let V-j-c denote the set of all positive semidefinite operators on a finite dimensional 
Hilbert space 7i with inner product (•, ■). Then, for any finite natural number n, any finite 
p > 0, and any finite n-tuple (Ai, A 2 , . . . , A n ) of elements of Vn, define 



The main result of this paper is the following: 

Theorem 1 For < p < 1, <E> P is a jointly concave function of its arguments. For p = 2, 
$ p is jointly convex. For p > 2, $ p is neither convex nor concave. 

We conjecture that $ p is jointly convex for 1 < p < 2. We state all of the theorems 
in a finite dimensional context, and some of our methods of proof explicitly involve this 
finite dimension. Nonetheless, the results themselves do not depend on the dimension, and 
therefore easily extend to the appropriate trace classes on an infinite dimensional Hilbert 
space. 

We note that the trace in Theorem 1 is essential; the asserted trace inequalities do 
not hold as operators inequalities. If they did, we would have, for example at p = 2 that 
(A 2 + B 2 ) 1 / 2 < A + B. This is of course not true in general for positive operators, as is 
well known and easily checked. 

We shall use Theorem 1 to derive a Minkowski type inequality for traces of operators 
on a product of three Hilbert spaces. To set this in perspective, recall that the Minkowski 
inequality says that for non-negative measurable functions / on the Cartesian product of 
two measure spaces (X, fx) and ( Y, v) , 



for p > 1, and that the opposite inequality holds for < p < 1. 

A direct analog of (1.2) holds for positive operators A on the tensor product of two 
Hilbert spaces Hi <8> Ti.2- To state it, let Tr\A denote the positive operator on H2 that is 
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given as a quadratic form by 

(v,TnAv) = ^(uj (g)v,A(uj ®v)) 
j 

where v G H2 and the Uj constitute an orthonormal basis of Hi- As is well known, the 
quadratic form on the left is independent of the choice of the orthonormal basis on the 
right. The operator Tr^A so defined is called the partial trace of A over Hi- It will be 
convenient, and generally clearer, in what follows to write Tr\ also to denote the usual 
trace on Tii for operators A on 7ii alone. 

The following is the tracial analog of (1.2): 

Theorem 2 Let A be a positive operator on the tensor product of two Hilbert spaces 
n 1 (g)H 2 . Then for all p > 1, 

(Tr 2 (T ri A) p ) 1/p < T ri ((Tr 2 A p ) 1/p ) (1.3) 

and inequality (1.3) reverses for < p < 1. 

Returning to (1.2), note that it has a trivial extension to functions of three (or more) 
variables. Though trivial, it has an interesting consequence. If one considers a non-negative 
measurable function f(x,y,z) on the Cartesian product of three measure spaces (X,fi), 
(Yu) and (Z, p), and simply holds z fixed as a parameter, one gets 

J Qf /Or, y, z)dvjd^ ^ < J J^J r(x, y, z)d^ '\v (1.4) 
pointwise in z for p > 1. Integrating in z then yields 

L Ux Uy V " Z)dU j dfi ) -j z jy{j x F(X ' V " Z)dfX ) /l>MP (L5) 

for p > 1, and of course the inequality reverses for < p < 1. 

Now, since (1.5) is an equality at p = 1, we get another inequality by differentiating 
(1.5) with respect to p at p = 1. This yields an entropy inequality. In fact, using the 
homogeneity of (1.5), we can normalize / so that it is a probability density. Recall that 
for any probability density p on any measure space (X, p), the entropy S(p) is defined as 

S(p) = - [ plnpdp. (1.6) 
Jx 
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We then denote various marginal densities of / as follows: 

f2,3(y,z)= f(x,y,z)dfj, fi,a(x,z)= f(x,y,z)dv h{z) = / / f(x,y,z)dfjdv 

JX JY JXJY 

Then the derivative of (1.5) at p = 1 is 

S(fi,s) + S(h, 3 ) > S(f w ) + S(h) (1.7) 

which is the strong subadditivity of the classical entropy; see [L75]. 

Now consider operators on the product of three Hilbert spaces, and a density matrix 
A; i.e., a positive operator on Hi ®H 2 ®H 3 with TrA = 1. The entropy S(A) of a density 
matrix A is defined by 

S(A) = -Tr(AlnA) . (1.8) 

The operator analog of (1.7) is the Lieb-Ruskai [LR] strong subadditivity inequality 
for the quantum mechanical entropy: 

5(^1,3) + S(A 2 , 3 ) > S(A h2 , 3 ) + S(A 3 ) 
where, in analogy with our notational conventions for marginal densities, we define 

^1,2,3 = A A 2}3 = TnA , A 3 = T n Tr 2 A 

and so forth. 

Thus, the differential form of Minkowski type inequality (1.7) is known to hold at p = 1 
for operators. It is therefore natural to enquire whether there exists an operator analog 
of the three-variable Minkowski inequality (1.7) for other values of p. Unfortunately, the 
methods at our disposal suffice to establish this only for < p < 1 and for p = 2. 

Theorem 3 Let A be a positive operator on the tensor product of three Hilbert spaces 
H!®H 2 ®H 3 . Then 

Tr 3 (Tr 2 (TnAr) 1/p < Tr h3 ((Tr 2 AP) 1/p ) (1.10) 



for p — 2 and, trivially, p = 1, while the reverse inequality holds for < p < 1. 
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This is, nonetheless, enough to imply the strong subadditivity (1.9): one simply takes 
the left derivative at p = 1. 

It is readily seen by considering block-diagonal matrices that the inequality of Theo- 
rem 3 implies the convexity of $ p for p = 2, and the concavity of <E> P for < p < 1. By the 
same token, (1.10) cannot hold in general for p > 2 since this would imply the convexity 
of <E>p for such p, and Theorem 1 precludes this. This is in contrast to Theorem 2, the 
Minkowski inequality for two spaces, which holds for all p > 1. 

The fact that there is such an easy passage from the Minkowski inequality in two 
variables to that in three variables may leave one surprised that there should be any 
difficulty in making the same passage with operators. But difficulty there is. In fact, even 
the simple version in Theorem 2 seems to require a more intricate proof than does the 
corresponding statement for integrals - which after all is simply the statement that the 
unit ball in LP is convex for p > 1. In fact, we know of no previous proof of Theorem 2. 

We emphasize that there is no operator analog of the pointwise inequality (1.4). That 
is, if we omit Tr 3 on both sides of (1.10), the result will be two operators on 7i 3 , and these 
two operators do not satisfy the corresponding operator inequality. 

We present a proof of Theorem 2 in Section II. Then in Section III we prove Theorem 

I. In Section IV, we recast Theorem 1 into an equivalent form, from which Theorem 3 is 
readily derived in Section IV. Section V contains a brief comment on a relation between 
the conjectured convexity for 1 < p < 2 and a very interesting trace inequality of Birman, 
Koplienko and Solomyak [BKS]. 

II. Proof of Theorem 2 

The following proof of Theorem 2 is given for matrices, but is easily extended to 
operators as the statement is dimension independent. 

Let A be a positive operator on Vux®Tiii ^ ne tensor product of two finite dimensional 
Hilbert spaces. Suppose first that p > 1. We proceed by duality. There is a positive 
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operator B in V H2 with (Tr 2 (B q ) 1 / q ) = 1 with 1/q + 1/p = 1 such that 

(Tr 2 (T ri A) p ) 1/p = Tr 2 (BTnA) = Tr lj2 ((I ® B)A) = 
y]{uj <S> Vj, (I ® B)Am ® Uj) = y^(ttj <8> -Bfj, ® Uj) 

i,3 i,j 

for any pair of orthonormal bases {ui} and {vj}. We now choose the {vj} to be a basis of 
eigenvectors of B, and let {Xj} be the corresponding eigenvalues. 
Then the right hand side above becomes 



^ Aj (ui <S> Vj, A{ui <8>Vj)) < 

i,3 

\ 1 /q / / \p\Vp 

3 i j ' ' 

^(^((u i ®v j ,A(u i ®v j )) > j ^ 



Next, by the spectral theorem, for each % and j, 

(ui <S> Vj,A(ui ® vj)) < ((ui <g> Vj, A p (ui <S> Vj))) 1 ^ 
Using this, one arrives at 

(Tr 2 (TnA) p ) 1/p <Y,(j2 (<«< ® v *> A *<& ® v i)>) ) = 
]T((^,Tr 2 ^» 1/p 

i 

Now we choose the {ui} to be a basis of eigenvectors of Tr 2 A p . Then 

Y J ((u i ,Tr 2 A p u i )) 1/p = J2(ui{Tr 2 A p ) 1/p Ui ) = Tr{Tr 2 A p ) 1/p 

i i 

and the desired inequality is proved for p > 1. Note that this part of the proof works for 
all p > 1, not only 1 < p < 2. 

Now suppose < p < 1, and define r = 1/p and £? = A p so that A = B r . Since r > 1, 
the inequality proved above says Tri[(Tr 2 B r Y^ r ) > (Tr 2 (Tri_B) r ) 1 ^ r . Rewriting this in 
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terms of A and p, and switching the roles of Hi and H2, one obtains the desired result for 
<p < 1. □ 

III. Proof of Theorem 1 

As before we give the proof for matrices. Consider first the case < p < 1. The proof 
in this case proceeds by reduction to a theorem of Epstein [E] concerning the function 



A^Tr({BA p B) 1 ' p ) 



on Vji where B is any given element of T>u- Epstein's theorem says that this function is 
concave for < p < 1. 

To apply this, consider first the case n = 2 in (1.1), and define 



A 



and 



Then 



But 



Now define 



a = 



Ai 
A 2 



/ 
/ 



A p + oA p a = 



A p + A p 
A p + A p 



A p + aA p a = 2 [ 'l±°) A'f^+^l^W 1 -" 



n± 



2 

I±a 



and observe that these are complementary orthogonal projections. Thus, 



2Tr((A p + A p ) 1/p ) = 2 1/p Tr((n + ^ p n + ) Vp ) + 2 1 / p Tr((U_A p U_) 1/p ) (3.1) 



Epstein's theorem, with A = A and B = U± now implies that each term on the right hand 
side of (3.1) is a concave function of A, which means that the left hand side is a jointly 
concave function of Ai and A 2 . This concludes the proof for n = 2. 
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One now easily iterates this procedure to obtain the result for all dyadic powers n = 2 fc , 
and hence for all n. 

To prove the convexity of $2 there are several way to proceed, but the simplest was 
pointed out to us by S. Sahi. Namely, let n be given and consider the block matrix A 
given by 



A 



A-l ... 
A 2 ... 



A n 



Then 

<f> 2 (A 1 ,A 2 ,...,A n )=Tr\A\ 

where \X\ is the usual operator absolute value; i.e., \/X*X. In other words $2(^-1, A 2 , . . . , A n ) 
is simply the trace norm of A, is therefore clearly jointly convex in Ai, A 2 , . . . , A n . 

Finally, we show that convexity fails to hold for p > 2. To see this, choose any pair 
Ai,A 2 E Vui and any vector v such that 

(v, {{A\ + A p 2 )/2)v) < (v, ((Ai + A 2 )/2fv) . (3.2) 

Note the strict inequality here. It is always possible to find such Ai,A 2 and v for p > 2 
since, for such p, X 1— > X p is not operator convex. 

Now let II„ denote the orthogonal projection onto the span of v, and let 11^- = I — H v 
denote its orthogonal complement. Then, for a large number A to be fixed below, put 

B = U V + XU^ . 

Then, if $ p were convex, we would have 

limsuppr p ($ p (t Al ^ A2 ,B^j - ^ p (tA u B) - ^<S> p (tA 2 ,B^J < . (3.3) 

However, for small t > 0, 



2 J V V 2 

TrB+^Tr^-r^^y^J + 0(t 2 *>) 
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and 




P \ 



Thus, 




Now taking A sufficiently large, this last term on the right is stricly positive by (3.2). This 
contradicts (3.3), and thus convexity does not hold - not even separately. □ 

IV. Corollary of Theorem 1 and Proof of Theorem 3 

A corollary of Theorem 1 is obtained by writing the partial trace as an average, and 
exploiting the convexity and concavity established above. Let A be a positive operator 
on Hi ®H2- Next, suppose the dimension of H2 is N, and fix some orthonormal basis 
{ei, e2, . . • , ejv}- With respect to this basis, define the self-adjoint unitary operators Uij 
and Vij on H2 by 



where the i and j are a distinct pair of indices. Let Q be the subgroup of the group of 
unitary operators on H2 that is generated by this family together with the identity. Each 
operator W in this group acts by 



for some permutation ir(-), and some map s : 1,2, ...,N 1— > 0,1. Thus, the size of the 
group is 2 N Nl, and the point about it is that any operator on H2 that commutes with 
every element of this group is necessarily a multiple of the identity on H.2- Then 



Uij — I Ei } i Ejj + Eij + Ej^ 



Vi = I- 2£ M 



Wej = (-l) s(j) 




(4.1) 



WeG 
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From here one easily 



Theorem 4 For p > 0, let the map *& P (A) from positive operators A on Hi ®'H 2 to IR + 
be given by 



y p (A) = T ri ((Tr 2 A*>) 1/p ) . 
Then this map is concave for < p < 1, convex for p = 2, and neither for p > 2. 



(4.2) 



Proof: We shall assume that the dimension of H 2 is iV so that we may apply the averaging 
formula introduced above. We then have 

Tn((Tr 2 ^) 1/p ) =N 1 '*- 1 Tr lt2 (j(JjTr 2 A*®I n2 > J ^ = 
Nl/P ~ lTr U(wAn J2( I ^ W *) AP ( I ^ W )) 1/P ) = 



n i/p-i( 1 



1 /p / / \ 1 /p 

Tr lj2 1 I {(I ® W*)A(I ® W)) p 



The result now follows directly from Theorem 1. □ 

Notice that the conclusion of Theorem 4 not only follows from Theorem 1, but also 
implies it. To see this, suppose that the A in Theorem 4 is block diagonal with 



A = 



A 1 ... 
A 2 ... 

... Ar, 



Then clearly 



^ P (A) = %{A u A 2 ,...,A n ) 



We remark that if our convexity conjecture turns out to be true for 1 < p < 2, then 
a proof along the same lines as the proof above of Theorem 4 will prove the conjectured 
tracial generalization of Minkowski's inequality. 
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Proof of Theorem 3 Suppose that the dimension of Hi is N. The left hand side of 
(1.10) can be written in terms of namely 

Tr 3 (Tr 2 (T ri A)P) 1/p = Tr lj3 (tt 2 (JjTnA <g> ^ = 

y p (JjT ri A®I H ^ 

where the pair of spaces in the definition of ^ p is taken to be H 2 and Hi <8> H 3 . 
Then by (4.1) and the convexity of ^ p established in Theorem 4, 

9 p (±T ri A<2>I Hl ) =^(J^ 52(I®W*)A(I®W)) < 

V / V • We g J 

' weg 



The last term above is 

i/p 

TV, „ ( Tr„ (7 7" » W*\ AP(T » ^ 



E rr 1)3 (Tr 2 ((/®W*)^(/®W))) 



weg 

_ , 1/P 

: — Tr, „ ( r 7" M7*^ (7V„ 7" 6?i ^ 

2 JV 



^ Tr i.3((J®W*)(rr2i4 p )(/®W0) ' = 
' weg ^ ' 

Hvn E T ri J(I®W*)(Tr 2 Av) 1/p (I®W)) = 

' user. V / 

Tr 3 T ri ((Tr 2 AP) 1/p ) , 



2 N 



which is the desired result. □ 

V. The BKS Inequality and the 1 < p < 2 Conjecture 

Birman, Koplienko and Solomyak [BKS] proved that for p > 1, and A and B positive 
semidefinite operators, 

Tr(B p - AP) 1 ^ > Tr(B-A) + (5.1) 

where X + denotes the positive part of a self adjoint operator X; i.e., X + = (X + \X\)/2. 
In (5.1), neither B nor A needs to be bounded, but it is assumed that [B p — A p )^ p is 
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trace class. Though the inequality in (5.1) is only one of several very interesting inqualities 
proved in [BKS], we refer to it here as the BKS inequality. 

The proof is in two parts, the first of which is to reduce consideration to the case 
B p > A p in which case one has B = (A p + C p ) 1 / p with C > 0. Then (5.1) becomes 

Tr{C + A) >Tr(A p + C p ) 1/p (5.2) 

for all A > and C > 0. It is (5.2) that interests us here. 
Clearly (5.2) can be rewritten as 

%(A,C)<%(A,0) + %(0,C) , (5.3) 

which is a subadditivity property of <E> P for all p > 1. Since $ p is homogeneous of degree 
1, subadditivity and convexity are the same thing. Thus for p = 2, (5.3) is a special case 
of the convexity of $ p proved in Theorem 1, and for 1 < p < 2, it would be a consequence 
of the conjectured convexity for these p. However, the BKS inequality holds for all p > 1, 
not only for 1 < p < 2. 

There is a simple proof of (5.2) for matrices. Let 



M± 

so that 



Ap/ 2 
±Cp/2 o 



Tr(M^M±) 1/p = Tr(A p + C p ) 1/p . 
On the other hand, the spectrum of M±M± is the same as the spectrum of M±M±, so 

Tr(A + C) =Tr(M±M^) 1/p . 

One computes 



A p ±J 
±J C p 



M±Ml = 

with J = A p l 2 C p l 2 . Since X ^ Tr{X x / p ) is concave for p > 1, one has that 



Tr(A + C)=Tr( 2 -J >Tr(A p + C p ) /P 
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A recent application of the BKS inequality, and a different proof of (5.1) that holds 
in the case of unbounded operators, can be found in [LSS]. 
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